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Abstract 

The behaviour of solutions to fourth order problems is studied through the decomposition 
into a system of second order ones, which leads to relaxed formulations with the intro- 
duction of measure terms. This allows to solve a shape optimization problem for a simply 
supported thin plate. 
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1. Introduction 

In this paper we study the asymptotic behaviour of solutions of fourth order elliptic 
problems on varying domains. 

This has been widely studied in the past, in the case of second order elliptic oper- 
ators (see for instance [6], [9], [7]). We will use such results decomposing fourth order 
differential equations into a system of second order ones. 

Given a bounded open set U in IR n , n > 2 and a function / e H _1 (n), the fourth 
order equation 

A 2 u = f in H _1 (fi) 

Au E Hj(fi) (1.1) 
u E Hj(fi) 

is linked to the model for the vertical displacement u of an thin plate, occupying a region 
U , simply supported on dU , subjected to a load / . Simply supported means that the 
boundary is fixed, but that the plate is free to rotate around the tangent to dU . For the 
general treatment of plate theory we refer to [15], [11], [4], [14], [5]. 

In particular we want to study the asymptotic behaviour of solutions when the 
domain varies. To this aim we will show (Proposition 3.1) that problem (1.1) is equivalent 
to the system of second order equations 

-Au = vin W\U) 
u E Hj(tf) 
-Av = f in H _1 (*7) 

venl(u). 

This problem can be handled with the theorems valid in the second order case (Theorems 
2.1 and 2.2), and it will be proved that if u n are the solutions of problems like (1.2) on 
a sequence of subdomains U n of a given bounded domain Q , then a subsequence of u n 
converges weakly in Hq(O) to a function solving 

-Au^ + \iu^ = v 

<veHi(Q)nL2(fi) 
-Av + [iv = f 
.GHj(O)nL^). 

where \x is a measure. 

The study will be carried on for general fourth order elliptic operators with constant 
coefficients and no lower order terms, that can be splitted into two second order ones. 



(1.2) 



(1.3) 
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A motivation for the study of the asymptotic behaviour of solutions of Dirichlet 
problems in varying domains without geometric assumptions on the domains U n are the 
so-called shape optimization problems: given a function j : O x IR —> IR we consider the 
following problem 

r ™n / j(x, Uu (x))dx, (1.4) 

U&A($l) J 

a 

where U{Q) is the family of all open subset of O and u v is the solution of the problem 
of type (1.1) in the set U . 

In section 6 it will be shown that, in general, problem (1.4) does not have a solution. 
Hence a relaxed optimization problem will be introduced, where the set over which we 
minimize is the set of functions , where \x is a measure and solves the relaxed 
problem (1.3). This set is the closure of {u v : U G U(0,)} in L 2 (0). This problem will 
always have solution and its minimum will coincide with the infimum of integral in (1.4). 

This shape optimization problem, in the second order case, was studied in [1], [2], 

[3]- 



2. Notations and preliminary results 

Given an open subset U of IR n , Hq([7) is the usual Sobolev Space, H _1 ({7) its dual, 
and (•, •) the duality pairing. If U C 0, and u G Hq([7) , then the function 



u := 



u in U 
in Q. \ U 



is in Hq(O). From now on we will always denote, with the same symbol «, a function 
and its extension u. 

In this paper we will deal with elliptic operators L : V'(Q) — > V'(Q) of the form 



Lu = ^ c a d a u, 



where a and (3 are multiindeces, and m is the order of the operator that will be 2 or 4. 
In any case they will be without lower order term and with constant coefficients. 
The operators will be assumed to be elliptic in the sense that 



E > Tier, veGM n , 
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where 7 is a real positive constant. This, in our case, is the same as 

P(0^o veeK n \{0}, (2.1) 

where P is the polynomial ^ c Q f Q associated to the operator L. 

a I =m 

In this work, differential problems are always meant to be solved in the usual weak 
sense. This means, for instance, that, for u G Hq(£7) the expression 

-Au = f in R-\U) 

is an equality of linear functionals 

(-Au,v) = ^2 diudivdx = (f,v), 
* u 

for any v G H^([7) . 

As said above, the limit of a sequence of solutions of Dirichlet problems is not, in 
general, the solution of a problem of the same kind, but is the solution of a problem 
where a measure term appears. To deal with these problems we need to recall some 
notions. 

For the notion of capacity of a set £cO, which we will indicate by cap (E) , we 
refer to textbooks as [12] or [13]. We shall always identify a function u G Hq(O) with its 
quasi-continuous representative. 

Now, let Mq{VL) be the set of Borel measures which are zero on the sets of zero 
capacity. 

For one such measure fi G A4q(Q) , L^([7) will be the space of functions such that 

\u\ 2 dfi < +00. 



u 



Given a second order operator A = ciijdidj , with constant coefficients, for a function 

ij 

u G Ho([/) H L^(U) to solve the equation 

Au + \iu = f, 

will mean that 

J ciijdjiidiV dx + J uv dfx = J fvdx, (2.2) 



2-1 u u u 
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for all test functions in v G Hq(U) D 1%(U) . 

It can be easily proved that the space Hq(?7) D Yr^U) is a Hilbert space whenever 
H is in jV(q(0) , and hence, by Lax-Milgram Lemma, we have existence and uniqueness 
of solutions for a problem of the form 



for any linear elliptic second order operator A. 

The decomposition of a fourth order problem in a system of two second order equa- 
tions allows us to study the asymptotic behaviour applying well known theorems for the 
second order case separately to each equation. The following results, that can be found, 
for instance in [9] , [7] , are the key points of the theory. 

Theorem 2.1. Let A be a second order elliptic operator, as described above. For every 
sequence {/U n } in M.o{Vt) there exists a subsequence (i nk such that, for every sequence 
{g n } in H _1 (0) , strongly convergent to g G H _1 (0) , we have 



respectively. 

Theorem 2.2. For any measure fx G A4q(Q) there exists a sequence U n of open subsets 
of Q such that, for any sequence {g n } in H _1 (0) , strongly convergent to g G H -1 (fi) , 
the solutions z n of the problems 




z nk — z weakly in H (Q), 



where z. 



and z solve 




Az n = g 
z n G Hj( 



'n 



in H" 1 



(Un) 



(Un) 



converge weakly in H, 



(O) to the solution z of 



( Az + fiz = g 
{«GHj(0)nLj(0). 
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3. Decomposition of fourth order operators 

Fourth order elliptic problems are studied mainly with two different kinds of bound- 
ary conditions. In the model case of the bi-laplacian, they correspond to two different 
physical problems regarding, as said above, the displacement of a thin plate. Problems 
of the type 

A 2 u = f in R-\n) 
u E Hg(fi) 

3u 

correspond to having u = — — = on <9fi , that is to say that the plate is clamped along 

on 

its boundary. The asymptotic behaviour of the solutions of such problems has been 
studied in [10]. 

In this work we deal with the second kind of boundary conditions, as in (1.1), and 
we do it decomposing the problem in a system of two second order equations. 
The decomposability of the fourth order operator 



Lu = c a d 



|a|=4 



can be seen through the associated polynomial 

\a\=4 

It is a simple algebraic fact that, if the polynomial can be splitted into two second 
degree polynomials 

then other decompositions can be obtained only by exchanging the order or multiplying 
and dividing by constants. Observe that, according to (2.1), if P is elliptic, then so are 
R and Q. 

We remark here that such a decomposition can always be done in the two dimensional 
case. In higher dimensions this is not always possible. 
So assume that 

n I n \ 

Lu = E c ^ a = E b v d 3 d i E a ki d k9i u = BAu. 

\a\=4 i,j=l \k,l=l I 
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Proposition 3.1. Let U be a subset of Q , and A = ^ ciijdjdi, B = / J bjjdjdj be 

second order elliptic operators with constant coefficients. Let f G H _1 ([7) . The following 
three problems are equivalent: 



( BAu = f in H" 1 ^) f U G H o (C/) : Au G L2(C/) (Au = v in H _1 (C7) 

(i) ^ e Hq([/) («) /^uB^dx = </,„> (in) U B ^fl H . 1(c0 

l ueH o^J lv^GH5(Z7) : B^L 2 (f/) [ i; G Hj(l7) 

(3.1) 

Proof, (i) =>- (ii) . The equality {BAu, <p) = (f, <p) holds in particular for any <p in 
B.q(U) such that G L 2 ([7). Since £? is symmetric and Au G Hg(i7) we get 

J AuBipdx = (Bcp,Au) = {BAu, ip) = (/, ip) 
u 

(ii) =>- (Hi) . Let v be solution of 

ven l (u), (Bv,<p) = (f,<p) W>gH (£/). 

Then 

(B<p,v) = {f,<p) 
and if, in particular, cp is such that G L 2 (J7) , we get 



J vBipdx = (/, (/?), 



which, subtracted to the equation in (ii) gives 

y (Au - v)Bcpdx = 0. 



u 



Observe now that, thanks to Lax-Milgram theorem, every function in L 2 ([7) can be 
written as Bip, with (p G Hq(U) , so we obtain 

(Au - v)z dx = 0, VzeL 2 ((7). 



?7 
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Taking z = Au — v, we get that \\Au — vW^^ = 0, hence Au = v, and u solves problem 
(Hi) . 

(Hi) =>- (ii) . If u is a solution of (in) then 

{Bv,<p) = {f i < P ) v v eKl(u). 

and if, in particular, Bip G L 2 (£/~) we have 

J vBipdx = (/, y>). 

Since now Am = v we have 

y AuBipdx = (/, </?), 

C7 

for every <p G Hq(C7) with G L 2 (L r ) , hence u solves (ii) . 

(ii) (i) . We have already proved that from (ii) it follows that Au G Hq([7), 
hence BAu G H _1 (t/") , and so 



(f,(p) = J AuBcpdx = (Bcp,Au) = (BAu,(p); 
u 

that is, BAu = f in H _1 (i7) □ 



Remark 3.2. The equivalence af these problems gives for free existence and uniqueness, 
since this is true for problem (Hi) thanks to the observation made in section 2. This 
is not obvious for problem (ii) because the spaces of solution and of test functions are 
different. 



Remark 3.3. Notice that if the boundary of U is regular, then by regularity theorems, 
Au G L 2 ([7) implies that u belongs to H 2 (L r ). The same is true for the test functions 
tp. So the problem(3.1) (ii) becomes: 

( u G Rl(U) f]R 2 (U) 
AuBip dx = (f, (p) 
Vv? G Hj(f/)nH 2 ((/). 
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Remark 3.4. It is important to remark that the boundary conditions in the equivalent 
problems (3.1) (i)(ii)(iii) depend on the choice of the decomposition L = BA. 
If we have two decompositions, in the sense that 

Lip = BxAxp = B 2 A 2 (p, Vv? e V(Q), 

then the boundary conditions in (3.1) (i) are different. Hence also the other problems 
(ii) and (Hi) differ. Had we chosen the Dirichlet boundary conditions, that is seeking 
u G Hq({7), all decompositions would have given the same solution. But in this case 
what makes the difference is the boundary condition, as can be seen with the following 
example. 

Example 3.5. The operator 

Lu . Uxxxx ^^"yyyy ^U X xyy 

can be written, for instance, as the product of 

Au = Au and Bu = Au + u yy . 

But the integral in equation (3.1) (ii) will be different according to which operator we 
will apply first: 

J AuBt) dx = J (u xx V xx -\- IllyyVyy -\- 2,Uxx^yy UyyV X x) ^X , 

j BUAV (IX = j (U XX V XX ~\- IllyyVyy ~\~ UxX^yy ~\~ 2lUyyV ' X x) dX . 

Computations show that these two integrals differ by a term on the boundary, which 
would vanish if functions where in Hq(LZ) . 

4. The asymptotic behaviour 

We come now to examine problem (3.1) (i) when we have a sequence of domains U n 
all contained in O. Let / be in H _1 (0) (this implies it is also in B.~ 1 (U n ) for any U n ) 
and consider always functions of B.Q(U n ) trivially extended to the whole of O. As proved 
in Proposition 3.1, we can study directly 

Au n = v n in H _1 ([/ n ) 

u n eRl(U n ) ( . v 

Bv n = / in Ur\U n ) 1 ' ) 

v n eR 1 Q (U n ). 
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We first consider the second equation 

Bv n = f inH" 1 ^) 
v n eK l (U n ). 

From Theorem 2.1 we know that there exist a subsequence, which we still call U n , and 
a measure /ig, depending on the sequence U n , on the operator B , but not on /, such 
that 

v n — ^ v weakly in Hq(O) 



and v solves 



r Bv + fi B v = f 

\^GHi(0)nLj B (0), ^ 



in the sense specified in (2.2). 

We can now apply Theorem 2.1 to the problem in u 

Au n = v n in H _1 (t/" n ) 

taking as U n only those in the subsequence obtained for B . Again there exists a subse- 
quence, which we still call U n , and a measure ha , depending on the sequence U n , on 
the operator A , but not on the sequence v n , such that 

u n — ^ u weakly in Hq(^) 

and u solves 

J Au + hau = v 

\«GHi(0)nL2 A (0). W 

This allows us to conclude that, if u n are the solutions of system (4.1) then, up to a 
subsequence, 

u n u weakly in Hq(O) 

and u is the solution of 

Au + hau = v 
^Hi(Q)nLj A (Q) 

Bv + hbv = / ' 
,GHj(0)nL2 fl (0). 
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5. The single equation formulation 

The goal of this section is to write problem (4.4) as a single equation of the form 

J (Au + u A u) (Bcp + u B ^>) dx={f,(p). 



n 



Of course Au + \iau alone doesn't make sense, because it has to be under- 
stood in the sense explained in chapter 2. What we will do now, is hence to define 
suitable function spaces, which will play the role of {z E Hq(O) : Az E L 2 (0)} and 
\w E Hq(O) : Bw E L 2 (0)} and operators A MA and A AtB which give meaning to the 
integral 

J A MA it A^ B ipdx. 
n 

The construction of such operators will be done only relatively to operator A , being the 
one relative to B perfectly analogous. 
Consider the operator 

R A iH^^Hj^nLjJf]) 

z i — ► w 

where w is the solution of 

j Aw + haw = z 

i^eHj(n)nLj A (n) ; C 5 - 1 ) 

so that we have u = Rav (and v = Rsf)- 
Define now the space 

V A (n) := {w ERl(Q) HLlJQ) : 3z E Hq(O) fl L 2 tA (0) s.t. w = R A z) . 

Let us prove that the function z is unique, for each w. 

Assume, by contradiction, that there exist two 21,236 Hq(Q) fl L 2 A (fi) such that 
w = Raz\ = Razi ■ From the equation we have 

J ( Zl - z 2 )^dx = 0, vv> e Hj(ft) n Lljn). 

Since z x - z 2 E Hj(O) n L 2 lA (^) we get 

Il 2 (h) 



/ 



(zi - z 2 ) 2 dx = \\zi - ^H^rn) = 
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and we conclude that z\ = zi . 
So we have defined 

v A (n) = R A (Hj(0)nLj A (0)) 

and on its image Ra is one-to-one. So, on the space Va(Q) , the inverse operator can be 
defined: 

A MC :F A (fi)-^Hj(Q)nl4(Q). 

w i — ► z s.t. w = Raz 
Observe now that the operator Ra is symmetric, that is: 

(h,R A g) = (g,RAh), 

for any h G H _1 (0) and g G H^(fi) . 

We are now able to write system (4.4) as a single equation. We have v = R B f 
and u = Rav, so that u G Va(O). For any ip G V B (Q) there exists a unique z G 
Hq(O) n L^ B (0) such that <£> = R B z. Hence, 



(f,<p) = (f 1 R B z) = (z,R B f) = (z,v) 

= (Amb^V") = J A t*B i P A^ A udx. 



Last equality holds because both A^ B p and A^ A u are in L 2 (Q). 

Hence we can conclude that the function -u, solution of (4.4), also solves 

u G V A (Q) 

J A MB ^> A^udx = (f,<p) 

yip g v B (n). 

It is easily seen that also the converse holds. Let v be a solution of 

Bv + \i B v = f 

^GHi(0)nL2 B (0) ; 

then, by the definition, v = R B f. Hence, taken any z G Hq(O) n L^ B (0), we set 
<p := R B z, and we have: 



(v,z)=(R B f,z) = (/, R B z) = (f,<p) 
= (A fIA u,A^ B p) = (A^ A u,z). 
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Taking then z = v — A MA w, we obtain 



so 



1^ - A MA«II L 2 =0, 



A MA « = v, 



that is, by the definition of A MA , 

u = Rav, 

hence u solves (4.4). 

All this can be summarized in the following 

Theorem 5.1. A function u G Hq(Q) flL^ A (0) is a solution of 

u G V A (Q) 

J A MB ^ A^ A udx = (f,(p) 
if and only if it solves 

Au + hau = v 
«GHi(0)nLj A (0) 
Bv + hbv = f 
,GHj(0)nL2 B (fi). 

It is importat to remark that the road back to the equivalence with a fourth order 
problem stops here. This is because, in general, the image of A MA is not contained in 
the domain of A^ B , even if A = B and the two measures coincide. 



6. The optimization problem 

Let us consider now the optimization problem (1.4). 

Let / be a function in H _1 (0) and let j : O x IR — > IR satisfy the standard 
Caratheodory conditions and be such that 

\j(x,s)\ < b(x) +P\s\ p , for a.e. x G O and Vs G IR, (6.1) 

2n 

with b G L x (0) , (3 G TR and 1 < p < 2* , where 2* = is the Sobolev exponent. 
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We want to study 

r ™n / j(x, Uu (x))dx, (6.2) 
ueu(Q) J 
a 

where W(0) is the family of all open subsets of and u v is the solution, trivially 
extended in \ U , of the problem 

A 2 u = f in H -1 (17) 
Au G Hl(U) 

In general problem (6.2) does not have solution, in the sense that the infimum is not 
attained on the set {u u : U G U(Q)} . This can be seen with the following example. 

Example 6.1. Consider the case of O C lR n with smooth boundary. Consider a function 

w G C°°(tt), s.t. w > in O, w = and Aw = on <9fi 

and set 



Consider now 



so that the integral 



z:=-Aw + w, f:=-Az + z. (6.3) 



j(x,s) = (s- w{x)f , 



2 

(u(x) — w(x)) dx (6.4) 



takes the value zero if and only if u = w a.e. It will turn out from the theory (see 
equation (6.7)) that zero is in fact the infimum of the functional over all solutionds of 
Dirichlet problems on subsets of O. 

We prove now that, with this choice of /, w can not be the solution of the problem 

A 2 u v = f in H _1 (17) 
Au v G HS(C0 
u v GHi(*7), 

for any U C O, or, equivalently (thanks to Proposition 3.1), of system 

-Au v = v in R-\U) 

-Av = f in R-\U) 
venl(U). 
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First consider the case when cap(0 \U) =0. Then the spaces Hq(O) and Hq(U) 
coincide (see for instance Theorem 4.5 in [13]), and the problem in U is the same as the 
problem in O. We show that w ^ uq. From (6.3) we have 

A 2 w - Aw = -Az 

= f-z 

= f + Aw — w. 

w 

So A 2 w = f + 2Aw — w, and hence w can be such that A 2 w = f only if -Aw = 

2 

but this is impossible because, by the maximum principle, w should be negative, while 
we have choosen it strictly positive. 

The second case is when cap(0\[7) > 0. If this happens u v has to be zero in Q\U , 
so that / (o on a set of non-zero capacity. But for functions in Hq(O) to be equal 
Lebesgue-a.e. is the same as cap-a.e., hence u v can not be a minimizer for (6.4). 

The example shows then, that to be able to solve always our optimization problem, 
it is convenient to seek our minima in the larger set M := {u^ : \x E Mo(fl)} , which is 
the closure of N := {u v : U E U(Ci,)} in the weak topology of Hq(O) . In other words 
we introduce the relaxed optimization problem 

min / j(x,Uu(x)) dx, (6.5) 

Q 

where w M denotes the solution fo the relaxed problem 

\ E « nl » (6.6) 
i; M GHS(fl)nLj(f2) 

and each equation is meant in the sense of formula (2.2). 

In order to see that solving the new problem (6.5) gives the complete solution to 
problem (6.2), we have to show two things: that (6.5) has a solution and that 

min / j(x, «m(x)) dx = inf / j(x,u u (x)) dx. (6.7) 
neM {Vt) J ueu(n)J 

As for the first thing, we need to show that the set M is closed in the weak topology 
of Hq (O) . So consider minimizing sequence [i n and let u^ n be the solutions of the 



Asymptotic behaviour of solutions of fourth order Dirichlet problems 



15 



systems 

-AVf, n + V, in fl n = f 

^eHj^nLjJfi). 

Applying Theorem 2.1 to the second equation we have that there exist a subsequence, 
called for simplicity fi n , and a measure fx G A4q(Q) such that v^ n — ^ v M weakly in Hq(Q) , 

-Av M + v^n = f 
VM GHi(0)nL2(0). 

Applying again Theorem 2.1 to the first equation (only those corresponding to the 
subsequence), we obtain that 

tVn w m weakly in Hq(O). 

and is the solution of (6.6). 

Hence is admissible and, by the minimizing property of the sequence {/U n }, 

/j(x,u fI (x)) dx = min / j(x, u v {x)) dx. 
.vi,,(<>s J 

n q 

On the other hand observe that every problem of the form 

— Az = w 
z e Rl(U) 
—Aw = g 
w e E&U) 

can be viewed as a relaxed problem with the measure 

U(r>\ /0 if B\U has capacity zero 

^ ' \ +oo otherwise 

since this holds for second order problems as can be seen, for instance, in [2]. 

This ensures us that problem (6.5) is indeed an extension of (6.2), in the sense that 



N C M , so that 



^} n ^ / j(x,Uf,(x))dx< inf / j(x,u u (x)) dx. 
.Mo(fi) J ueu{ci) J 



a, n 
To prove the inverse inequality we need the following 
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Lemma 6.2. Let \x G A4q(Q) . Then there exists a sequence U n of subsets of O such 
that 

u Un — ^ weakly in Hg(O) 

where u Un and solve respectively 

-A Uu = v n in B.-\U n ) ( -Aw M + u^fi = 

\ eM( ^^rr, anJ\ eH ° (!i)n H (!1) 
-Av n = / m H \U n ) j -Au M + v^ = f 



Proof. We know, from Theorem 2.2, that given we have a sequence of sets t/ n such 
that, for every g G H _1 (0) , if 

-Az n = gin W\U n ) 
z n G Hj(t/ n ). 

then the sequence z n tends weakly in Hg(O) to z M , solution of 

z M GHS(n)nLj(n). 

Since the sequence U n depends only on the operator and not on the right hand side the 
same sequence of sets suits for the problem in u . Hence u v converge to w M solution of 
the second system. □ 
Since, by the hypothesis (6.1), the operator u » J j{x ,u)dx is continuous in the 

strong topology of L 2 (0) (see [2] Theorem 4.1), the density result enables us to conclude 
that 

min / j(x,Un(x)) dx = inf / j(x, ua(x)) dx. 
veM (n) J u&a{u) J 

This ensures us also that every solution of a relaxed system can be approximated by 
solutions of fourth order problems. 

Remark 6.3. We remark here that the model case of the bi-laplacian can be extended 
without changes in the proofs to the case of a fourth order elliptic operator with constat 
coefficients L such that 

Lu = C 2 u, 

where C is a second order elliptic operator with constant coefficients. 
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